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A partial wave analysis (PWA) of tlie of the r;7r° system (where r) 77) produced in the charge 
exchange reaction vr^p rjTv^n at an incident momentum of 18 GeV/c is presented as a function 
of rjiv" invariant mass, m^^o, and momentum transfer squared, from the incident 7r~ to 

the outgoing jyTr" system. S, P and D waves were included in the PWA. The ao(980) and 02 (1320) 
states are clearly observed in the overall Tyvr" effective mass distribution as well as in the amplitudes 
associated with S wave and D waves respectively after partial wave decomposition. The observed 
distributions in moments (averages of spherical harmonics) were compared to the results from the 
PWA and the two are consistent. The distribution in t^-^^^ for individual D waves associated with 
natural and unnatural parity exchange in the t-channel are consistent with Regge phenomenology. 
Of particular interest in this study is the P wave since this leads to an exotic J^'^ = 1 for the 
r]-n system. A P wave is present in the data, however attempts to describe the mass dependence 
of the amplitude and phase motion with respect to the D wave as a Breit-Wigner resonance are 
problematic. This has implications regarding the existence of a reported exotic J^'~^ — meson 
decaying into riTv^' with a mass near 1.4 GeV/c^. 

PACS numbers: 13.20. Jx;14.40.Cs 
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I. INTRODUCTION 

The observation nearly four decades ago that mesons 
are grouped in nonets, each characterized by unique val- 
ues of J^*^ - spin (J), parity (P) and charge conjugation 
(C) quantum numbers - led to the development of the 
quark model. Within this picture, mesons are bound 
states of a quark (g) and antiquark (g). The three light- 
quark flavors {up, down and strange) sufRce to explain 
the spectroscopy of most - but not all - of the lighter- 
mass mesons (i.e. below 3 GeV/c^) that do not explicitly 
carry heavy flavors (charm or beauty). Early observa- 
tions yielded only those J^^ quantum numbers consis- 
tent with a fermion-antifermion bound state. The J^'~" 
quantum numbers of a qq system with total quark spin, 
S, and relative angular momentum , L, are determined 
as follows: J^L + S,P= (-1)^+^ and C = (-1)^+^. 
Thus J^^ quantum numbers such as , 0^ , 1 and 
2^ are not allowed and are called exotic in this context. 

Our understanding of how quarks form mesons has 
evolved within quantum chromodynamics (QCD) and we 
now expect a richer spectrum of mesons that takes into 
account not only the quark degrees of freedom but also 



the gluonic degrees of freedom. Gluonic mesons with no 
quarks {gluehalls) are expected. These are bound states 
of gluons and since the quantum numbers of low-lying 
giueballs are not necessarily exotic, they should manifest 
themselves as extraneous states that cannot be accom- 
modated within qq nonets. Indeed there is evidence for 
a glueball state 1] based on overpopulation of the scalar 
nonet but identification is complicated since the observed 
states can be a mixture of a glueball and qq. Excitations 
of the gluonic field binding the quarks can also give rise 
to so-called hybrid mesons that can be viewed as bound 
states of a quark, anti-quark and valence gluon (qqg). An 
alternative picture of hybrid mesons, one supported by 
lattice QCD, is one in which a gluonic flux tube forms 
between the quark and anti-quark and the excitations of 
this flux tube lead to so-called hybrid mesons 2]. Con- 
ventional qq mesons arise when the flux tube is in its 
ground state. Some hybrid mesons can have a unique sig- 
nature, exotic J^'" , and therefore the spectroscopy of ex- 
otic hybrid mesons is not complicated by possible mixing 
with conventional qq states. According to the flux-tube 
model 01 and lattice gauge calculations Q, one expects 
the lightest J^'^ = 1 exotic hybrid to have a mass of 
about 1.9 GeV/c^. In contrast, calculations based on the 
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MIT bag model ^ |H IS place the mass of the lightest 
J^'-^ = 1 ^ hybrid meson at about 1.4 GeV/c^. 

The spectroscopy of exotic mesons provides an attrac- 
tive starting point for the study of gluonic excitations. 
Since flux tubes are thought to be responsible for the 
confinement mechanism in QCD, experimental informa- 
tion on this spectroscopy is of fundamental importance. 



A. Experimental situation 

After about two decades of experimental searches there 
have been reports of experimental observations of states 
with J^^ — 1 ^ by the Brookhaven E852 collabora- 
tion in TT^p interactions at 18 GeV/c. One of these has 
a mass of (1593 ± St^S) MeV/c^ and width of (168 ± 
20tll°) MeV/c^ and decays into p°7r" HQ and another 
has a similar mass, (1597 ± lOt^o) MeV/c^, but a larger 
width,(340 ± 40;^|^) MeV/c^, and decays into ry'Tr" 0. 
The E852 collaboration also reported observation of an- 
other J^^ = 1-+ state with mass (1370±16t^°) MeV/c^ 

and a width of (385±40j;io5) MeV/c^ decaying into rj-K' 
[lOl ITU . If an rjn system is in a P wave, the result- 
ing J^'~^ quantum number combination is exotic (1 ^). 
Critical to the identification of this state is not only 
showing the presence of a P wave, but also that the re- 
sulting line shape is consistent with a Breit-Wigner and 
that the phase motion of the P, as determined by its 
interference with the dominant D wave cannot be due 
solely to the a^(1320) resonance. Soon after the E852 
report, the Crystal Barrel Collaboration reported an ex- 
otic J^'-^' = 1 ^ state produced in pn — > tt^tt^t; obtained 
by stopping antiprotons in liquid deuterium (12] . They 
reported a mass of (1400 ± 20j:^^) MeV/c^ and a width 
of (310±50j:^JJ) MeV/c^. A later analysis of^ Tr°n°r] 
supported this evidence for an exotic state . 

The first claim of an exotic meson decaying into rjir^ 
with a mass of 1400 MeV/c^ was made by the GAMS 
collaboration in the reaction n^p — + rjir^n |14| but a later 
analysis by the group ^o ambiguous results. The 

VES collaboration also presented evidence for a P-wave 
contribution in r/7r (16j] and at KEK a claim was made 
for an exotic rjn state |l7j as well but with a mass and 
width close to that of the 02 (1320) and leakage from the 
dominant D wave could not be excluded. 

Indeed, in all the observations, the rjn P wave enhance- 
ments have cross sections that are substantially smaller 
than the dominant 02 (1320) so this leakage, usually due 
to an imperfect understanding of experimental accep- 
tance, is a source of concern. For example, an observed 
enhancement in the P wave in the pir system at low mass 
was reported in references 0,13 as being due to leakage. 
Apart from these experimental issues, the interpretation 
of the nature of low-mass rjir P wave amplitude and phase 
motion should be guided by the principle of parsimony - 
less exotic interpretations must also be considered. 

This work presents a partial wave analysis (PWA) of 



the rjTr'^ system produced in n^p interactions at an inci- 
dent momentum of 18 GeV/c using the E852 apparatus. 
In this analysis both the 77 and 7r° are detected through 
their 27 decays. The starting point uses data correspond- 
ing to the reaction tt^p Ajn. A subset of these data 
correspond to n^p —> Tr'^n'^n and the partial wave analy- 
sis of the TT^TT*^ system has been reported on in an earlier 
publication 0| . Here we present the PWA of data from 
TT~p — > jyTT^n. 

In contrast to the r]TT^ system, charge-conjugation C 
is a good quantum number for the neutral rjTr^ system. 
In the E852 experiment the r]Tr~ and rjir'^ systems have 
different systematics. The latter relies on measurements 
from an electromagnetic detector alone while the former 
also requires information from charged particle tracking 
in the forward direction. Moreover, as will be discussed 
below, two resonances dominate the 7r~p — > rjir'^n reac- 
tion; the scalar {S wave) ao(980) and the tensor (D wave) 
02(1320) providing benchmarks for the PWA. In con- 
strast, only the 02 (1320) is dominant in the Tr~p — > rin~p 
reaction. 



B. Outline for this paper 

This paper is organized as follows. The experimental 
overview is presented in Section 2. Event reconstruction 
and data selection are described in Section 3. The details 
of the PWA formalism and results are given in Section 
4 along with the criteria used to select the physical so- 
lution from the set of ambiguous solutions. In Section 5 
the results of leakage studies are presented. In Section 6 
the distributions in moments (averages of spherical har- 
monics) are compared with moments calculated from the 
PWA solutions. Fits to the moment distributions are also 
presented. In Section 7 the mass and t dependence of the 
PWA solutions are presented and discussed. In particular 
the t distributions for the D waves are described within 
the context of Regge phenomenology. The conclusions 
are summarized in Section 8. 



II. EXPERIMENTAL OVERVIEW 

The E852 apparatus ^3 '^^s built around and in- 
cluded the Multi-Particle Spectrometer (MPS) at BNL. 
The data used for the analysis reported in this paper 
were collected in 1994 and 1995 using a nearly pure (> 
95%) beam of negative pions of momentum 18.3 GeV/c. 
A 30-cm-long liquid hydro gen target was surrounded by 
a cylindrical drift chamber |20| and an array of thallium- 
doped Csl crystals 21| arranged in a barrel, all located 
inside the MPS dipole magnet. Drift chambers were used 
to track charged particles downstream of the target. Two 
proportional wire chambers (PWC's), downstream of the 
target, were used in requiring specific charged particle 
multiplicities in the event trigger. A 3000-element lead 
glass detector (LGD) ^22j measured the energies and posi- 
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tions of photons in the forward direction. The dimensions 
of the LGD matched the downstream aperture of the 
MPS magnet. Photons missing the LGD were detected 
by the Csl array or by a lead/scintiUator sandwich array 
(DEA) downstream of the target and arranged around 
an aperture to allow for the passage of charged particles. 

The first level trigger required that the unscattered 
or elastically scattered beam not enter an arrangement 
of two small beam- veto scintillation counters located in 
front of the LGD. The next level of trigger required that 
there be no signal in the DEA and no charged particles 
recorded in the cylindrical drift chamber surrounding the 
target or in the PWC's (an all-neutral trigger). In the 
1994 run all layers of the cylindrical drift chamber were 
used in the trigger requirement whereas in the 1995 run, 
only the outer layer was used. A common off-line analysis 
criterion required no hits in the cylindrical drift chamber. 
The final trigger requirement was a minimum deposition 
of electromagnetic energy in the LGD j22| corresponding 
to 12 GeV. 

The LGD is central to this analysis. The LGD was ini- 
tially calibrated by moving each module into a monoener- 
getic electron beam. Further calibration was performed 
by adjusting the calibration constant for each module 
until the width of the tt^ and i] peaks in the 77 effective 
mass distribution was minimized. The calibration con- 
stants were also used for a trigger processor that did a 
digital calculation of energy deposited in the LGD and 
the effective mass of photons striking the LGD 22] . A 
laser-based monitoring system allowed for tracking the 
gains of individual modules. 

Studies were made of various algorithms for finding 
clusters of energies deposited by photons, including is- 
sues of photon-to-photon separation and position finding 
resolution. These are also described in reference 1231. 



III. EVENT RECONSTRUCTION AND DATA 
SELECTION 

Experiment E852 took data in 1994 and 1995 with an 
approximate seven-fold increase in statistics collected in 
1995 compared to 1994. 

The combined data sets taken in 1994 and 1995 contain 
approximately 70 million all- neutral triggered events. Of 
these events, approximately 13 million were found to have 
four photons in the LGD. Figure ^ shows a scatterplot of 
one di-photon effective mass against the other di-photon 
combination in the event. There are clear indications 
of the presence of tt^tt'^ and rjir^ events in the sample. 
The single clustering at =0.135 GeV/c^ is 

dominated by tt'^tt'^ events and the two clusterings at m. 
=0.135 GeV/c2 
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=0.540 GeV/c are dominated by 
rjTT^ events. 

Studies were carried out comparing the 1994 and 1995 
data sets for both tt^tt^ and rjir'^ events. In each case 
distributions in the decay angles of the tt'^tt^ (or ■qir'^) 
system were studied in mass bins of 40 MeV/c^ and in 
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FIG. 1: The plot of pairs of di-photon effective masses {mij 
vs. rriki) for all pairs of photons k, I) is dominated by the 
si gnal. Clear evidence is also seen for the production of 



four bins of t, where t is the square of the momentum 
transferred from the beam tt^ to the outgoing 47 system. 
A Kolmogorov test measured the probability that the 
two distributions (1994 and 1995 data) in each case were 
consistent with having been drawn from the same parent 
distribution. The conclusion of this study is that the 1994 
and 1995 data sets are statistically indistinguishable. 

Events consistent with the production of two 7r° 's dom- 
inate the sample. The partial wave analysis of the tt'^tt'^ 
events is presented in an unpublished thesis of J. Gunter 
[23j | and in a previous paper The identification of 

events from the reaction Tr~p — > TyTr^n is described in the 
unpublished thesis of R. Lindenbusch . 

The sample of 45,000 n^p rjir^n events was selected 
from the 13 million four photon events by imposing vari- 
ous analysis criteria. It was required that no charged par- 
ticles were registered in the MPS drift chambers or the 
cylindrical drift chamber surrounding the liquid hydrogen 
target. This cut was performed at the pre- reconstruction 
level and was based on total hit multiplicity. Any event 
with a photon within 8 cm of the center of the beam hole 
or the outer edge of the LGD was removed. 

Monte Carlo studies were carried out to determine the 
criteria needed to select r^ir^ events and eliminate 
backgrounds from other processes, especially from tt'^tt'' 
events. Events corresponding to tt'^tt", ryTr", 7777 and 
were generated and passed through detector simulation, 
reconstruction, kinematic fitting and other event selec- 
tion software. The returned from kinematic fitting to 
the 7rp rjTr'^n reaction hypothesis was required to be 
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FIG. 2: Distribution in the (a) r]Tr° effective mass; (b) 
momentum-transfer-squared ( ~t^-^^^) from the incoming 
TT^ to the outgoing rjTv'^ system. These distributions are not 
corrected for losses due to acceptance. 



less than 7.8 (95% C.L. for a three-constraint fit). The 
two-constraint for the tt^tt'^ hypothesis (no require- 
ment that the missing mass be consistent with the neu- 
tron mass) was required to be greater than 100. It was 
found that this latter cut eliminated true tt'^tt'^ events 
that have a poor fit to the irp tt^tt^ti hypothesis. A 
further demand was that none of the other final state hy- 
potheses considered {rjrin, ryVn) had a better than 
that used to select the rjir'^n hypothesis. 

The final criterion was that the Csl detector registered 
less than 20 MeV, a cut which ehminated events with 
low-energy tt'^'s. Background studies estimated the 
non-77 background in the np — s- ■qir'^n. The signal-to-noise 
ratio varies from about 5.1 : 1 in the t/tt" mass region 1.2 
to 1.5 GeV/c^ to 2.5 : 1 in the 7777" mass regions below 
this range and above this range up to 1.8 GeV/c^. 

The 777r'^ effective mass distribution is shown in figure 
|2Ia). The ao(980) and a2(1320) states, one a scalar and 
the other a tensor, are clearly observed. The presence 
of both states is important for this analysis since it pro- 
vides a basis, as will be discussed below, for selecting the 
physical solutions from among mathematically ambigu- 
ous solutions in doing the partial wave analysis. This 
distribution in riir'^ effective mass is not corrected for the 
acceptance of the apparatus but the PWA procedure does 
take into account the effects of experimental acceptance. 

The distribution in t^-^^j^ is shown in figure [^^b). 
Figure |31 shows the four-photon effective mass distribu- 
tion for events with a di-photon mass combination con- 
sistent with a TT*^ and the other di-photon mass either 
below or above the r/ mass region for a portion of the 
total four-photon sample. The resulting 777r° spectra do 
not show the enhancements at the ao(980) and 02(1320) 
observed in the t^tt" spectrum. 

A possible source of background in the ?77r'' in the 
ao(980) mass region are events corresponding to produc- 
tion of the final state ■qTr^Tr'^n and in particular that sub- 
set corresponding to /i(1285)n where the /i decays to 
oott''. There events can enter the sample under study if 
one of the tt^'s escapes detection. From a Monte Carlo 
study of this sample based on a measurement of 
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FIG. 3: Distribution in the four photon effective mass for 
events with a reconstructed 7r° and the other 77 effective 
mass combination falling below (solid) and above (dashed) 
the r] mass region. Neither distribution shows the structure 
observed for rjTv" events. 
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FIG. 4: Acceptance as a function of (a) jyvr" effective mass; (b) 
momentum-transfer-squared ( —t^-^^^) from the incoming 
7r~ to the outgoing t^tt" system. 



/i(1285) production in this experiment [T^ it is esti- 
mated that at most 10% of the events in the ao(980) 
region are due to this background process. 

The acceptance in 7777° effective mass and t^-_,^^ was 
estimated by generating Monte Carlo events correspond- 
ing to the reaction ir^p — > rjir'^n. The acceptance func- 
tions as a function of rjir'^ effective mass are shown in 
figure^ The acceptance in tjr-^jyTr increases for —t near 
zero. This comes about because of the all-neutral re- 
quirement in the trigger. The recoil neutron produced in 
the reaction can interact and cause a signal in the Csl 
detector thereby vetoing true events but for small 

values of \t\ the neutron cannot escape the LH2 target 
and the event will not be vetoed. 
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FIG. 5: Acceptance in Gottfried- Jackson angles for the jyvr" 
effective mass region from 1.3 GeV/c^ to 1.4 GeV/c^ for: (a) 
cos{6g.j)- (b) (j)GJ- 



IV. PARTIAL WAVE ANALYSIS 

Partial wave analysis is used to extract production am- 
plitudes (partial waves) from the observed decay angu- 
lar distributions of the rjir^ system. A process such as 
'K~p — > riTT^n, dominated by i-channel meson exchange, 
is conveniently analyzed in the Gottfried- Jackson refer- 
ence frame. The Gottfried- Jackson frame is defined as a 
right-handed coordinate system in the center of mass of 
the produced rjir^ system with the z — axis defined by the 
beam particle momentum and the y — axis perpendicu- 
lar to the plane defined by the beam and recoil neutron 
momenta. The decay angles (0, 4>) are determined by the 
r] four- momentum vector. At a fixed beam momentum, 
an event is fully specified by (m^,r, ^, 0)- The accep- 
tances of the E852 apparatus for the Gottfried- Jackson 
angles were studied using Monte Carlo simulations. Re- 
sults for the TyTr" effective mass region from 1.3 GeV/c^ 
to 1.4 GeV/c^ are shown in figure |5| 

The production amplitudes and their relative phases 
are extracted from the accumulated angular distribu- 
tions, for bins in m,,7r and t, using an extended maximum 
likeUhood fit to the distributions in {0,(j)) 23, 24., ,25j . 
Fits were performed for three bins in (or sim- 

ply t) which we label as \ow-\t\ for \t\ < 0.14 (GeV/c)^, 
medium-|t| for 0.14 < \t\ < 0.31 (GeV/c)^ and high-|i| 
for \t\ > 0.31 {GeV/cf. For each t bin the PWA fits were 
performed in 0.04 GeV/c^-wide t^tt" mass bins. 

Events corresponding to the reaction w^p —> rjir^n 
were processed through PWA software that extracts par- 
tial wave amplitudes and phase differences between var- 
ious partial waves. This procedure takes into account 
information about the acceptance of the apparatus. The 
PWA formalism used in this analysis is model indepen- 
dent since the three-body final state follows from a quasi- 
two body interaction. This is in contrast to the isobar 
model needed for more complex final states, for example, 
in which one makes specific assumptions about a series 
of sequential decays leading to the final state under con- 
sideration. 

The decay angular distributions are described by an 



expansion in partial waves that are products of a pro- 
duction amplitude for a state of angular momentum L 
and associated magnetic quantum number M and a de- 
cay amplitude, also for a given L and M . The angular 
momentum L is between the rj and the 7r°. The pro- 
duction amplitude takes into account the spin degrees of 
freedom of the initial and final nucleons and the natu- 
rality of the exchange particle. That is, the assumption 
is made that in this peripheral process the produced rjir'^ 
state results from a ^-channel exchange from the nucleon 
vertex to the meson production vertex. Alternatively one 
can view the production vertex as a time-reversed decay 
of the produced "qir^ state into the beam particle (tt) and 
the exchanged particle. 

Parity conservation in the strong interactions implies a 
relationship between the rnagnetic quantum number M 
and naturality. Following [23 we choose a basis charac- 
terized by the naturality of the exchanged particle. Natu- 
ral parity exchange obtains if parity and spin are related 
as P = (—1)'^ (e.g. p meson exchange) and the exchange 
is unnatural \i P = (—1)'^+^ (e.g. tt meson exchange). 

This basis, the reflectivity basis, is one that does 
not possess the mixed symmetries of the D-functions or 
spherical harmonics under rotations about the y axis. For 
M > the linear combination with positive reflectivity, 
corresponding to natural parity exchange, is: 

+i?^,o(0,^,O) ^ i=[i?i,o(0,e,O)-(-l)^i?^Mo('/',^,O)] 

(1) 

and the linear combination with negative reflectivity, cor- 
responding to unnatural parity exchange, is: 

-Dl,^{cj,,6,Q) EE -^[Di,,i<l>,9,0)+{-l)^'D^^,,{c^,e,0)] 

(2) 

For M = the reflectivity is negative and: 



-D^o(0, 0,0) EE D^o (0,0,0) 



(3) 



by 



The distribution of events in m^^, t, 9 and </> is given 



/(m^,,t,0,0) - |^A+^,,(m^,,<)+i^t/o(^,'/'>0)|2 

LM 

+ lE^LMK.,t) -i?^;o(0, 0,0)1^(4) 



LM 



where the A^m are the production amplitudes and we 
refer to the various waves as 5*0 = Aqq, P+ = Al^^ _P_ = 



^10' 



D+ = A^-y, D_ = A21, and Dq = A. 



20- 



In this analysis we considered L = 0, 1, 2 (the nomen- 
clature is 5*, P and D respectively). Waves are further 
characterized by a subscript if M = and -I- or — if 
\M\ = 1. As noted in 25], amplitudes with \M\ > 1 
are not expected to be important for the production of 
mesons in irp quasi-two-body processes. This assumption 
is borne out by PWA fits to these data indicating little 



6 



(a) 



oil! 



> 

u 
O 
o 

3 5000 
p 



iiiiij 



(b) 



1.2 1,4 1.6 1.8 



1.2 1.4 1.6 l.S 



T|7C eifective mass 

FIG. 6: PWA solution for the (a) So and (b) D+ waves as a 
function of jjtt" effective mass. All mathematically ambiguous 
solutions are shown. 



contribution from \M\ > 1. We present these results af- 
ter the discussion of fitting the moments distributions. 
The amplitudes considered in this analysis, Sq, Pq, P_, 
Dq, and D_ are produced via unnatural exchange and 
P_i_ and are produced via natural exchange. Waves 
of different naturality do not interfere. Since the nucleon 
polarization is not being measured, in equation Q the 
dependence on the nucleon helicity is implicit in the pro- 
duction amplitude, e.g. jD+p = jDj+P + I^J"!^ with 
the upper labels referring to the target and recoil nucleon 
helicities and the two amplitudes representing nucleon 
helicity flip and non-flip. 



A. Ambiguities and finding the physical solution 

There are multiple discrete sets of partial wave am- 
plitudes that can give rise to exactly the same angular 
distribution psj . It can be shown that for the production 
of two pseudoscalars with only S, P and D waves there 
are up to eight ambiguous solutions. For a given bin in 
m^^o these multiple solutions can be found by randomly 
seeding the fit (using sets of different initial values for 
the amplitudes) and doing multiple fits. 

It is also possible to find the ambiguous solutions an- 
alytically, starting with one of the solutions found from 
the fitting procedure using the method of Barrelet zeroes 
[2^ . In this analysis the ambiguous solutions found from 
multiple randomly seeded fits were found to be consistent 
with the ambiguous solutions found analytically. 

As an example, all ambiguous solutions are shown for 
the So and Dj^. waves for the full range in t in figure 
El The presence of the dominant 02 (1320) in the D 
wave leads to a clustering of solutions about the expected 
Breit-Wigner lineshape. The solutions for the S wave are 
more widely distributed. 

When possible, mass-bin to mass-bin continuity along 
with other independent criteria must be used to select 
the physical solutions among the ambiguous solutions. In 
this analysis the presence of the well-established qq (980) 
and 02(1320) states, as observed in the rjn mass distri- 
bution, provides a means by which to select the physical 



solution. For each m^^r mass bin the solution selected as 
physical is one where the Sq and 13+ amplitudes and the 
phase difference, A^{Do — Sq), all as a function of 7777 
mass, are simultaneously most consistent with two inter- 
fering Breit-Wigner amplitudes with masses and widths 
for the ao(980) and 02(1320) as listed in the Review of 
Particle Properties l27|- For example, referring to fig- 
ure El the Sq amplitude of the physical solution in the 
1.0 GeV/c^ region should be large and its corresponding 
amplitude small and in the 1.3 GeV/c^ mass region 
the amplitude should be large and the So amplitude 
should be small. It was found that the requirement that 
the So and solutions simultaneously lie close to the 
expected line shapes leaves a few ambiguities that are 
removed when the criterion of proximity to the expected 
value of A$(Z3o — Sq) is imposed. We note that although 
the 02(1320) is dominant in the D+ wave, the criterion 
for selection involved the phase difference between the Dq 
and 5*0 waves - the D+ and So waves do not interfere. In 
summary, these criteria were sufficient to uniquely select 
the physical solution. 

In the discussion below that follows the description of 
fits to experimentally observed moments distributions, 
we discuss additional quantitative measures for discard- 
ing other ambiguous solutions as well as the effect of their 
inclusion on the stability of the final results. In addition, 
as discussed below, the ^-dependences of the selected D- 
wave solutions are well-described by Regge phenomenol- 
ogy as is the ratio {Dq — D-)/D+, giving us confidence 
in the selection procedure. 

In figure [7| the selected physical solution in each of the 
low-|i|, medium- |i I and high-|t| regions is shown as a func- 
tion of rjir'^ mass for the So-wave, the wave and the 
A$(S'o — Do) phase difference. Also shown are the cor- 
responding line shapes and phase differences assuming 
two interfering Breit-Wigner forms for the oo(980) and 
02(1320) - the basis for the criteria used to select the 
physical solutions. We note that the intensity and phase 
motion of the S'o-wave is also consistent with the presence 
of both the oo(980) and the ao(1450). The oo(1450) res- 
onance parameters have been extracted from a coupled 
channel analysis of final states producted in pp annihi- 
lations [25. This may account for the excess of events 
in the 1.4 GeV/c^ mass region for the S'o-wave intensity. 
However because of the large statistical errors in this re- 
gion we have not been able to unambigously extract the 
ao(1450) resonance parameters. We will return to this 
issue in Sections VI and VII. 



V. LEAKAGE STUDIES 

Incomplete or incorrect understanding of the experi- 
mental acceptance can lead to serious systematic effects 
in this type of analysis. In particular, it is possible that 
part of a large, dominant partial wave (for example !?-(-) 
could appear as a spurious signal in a different partial 
wave. We refer to this phenomenon as leakage and study 
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FIG. 7: Physical PWA solution for the (top) So (middle) D+ 
waves and (bottom) A<l>(5'o — Do) phase difference for low-|f|, 
medium-|t| and high-|f| as a function of r^Tr" effective mass. 
The curves are the Breit-Wigner line shapes used as criteria 
for selecting the physical solution and more details are given 
in the text. 
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FIG. 8: Results of leakage studies. A sample of Monte 
Carlo events was generated corresponding to production 
of the 02(1320) populating D waves. The sample was 
passed through cuts simulating detector acceptance and then 
through the PWA fitting software. In one case (a) the detec- 
tor acceptance was used in the fitting and in the other (b) the 
acceptance information was ignored. 



A$(P+ — D+) phase difference is a strong function of rjir'^ 



it by intentionally using an incorrect model of the accep- 
tance. 

A Monte Carlo sample of pure 02(1320) events was 
generated including both Dq and D+ waves. This Monte 
Carlo sample included the effects of experimental accep- 
tance. This sample was analyzed in two ways, in one 
case our best estimate of the acceptance was used, in the 
second no correction for the acceptance was made. 

The analysis using the acceptance correction should 
produce as output only Dq and waves with all oth- 
ers consistent with zero. This was found to be the case 
and the resulting So and P+ solutions are compared to 
the intensity in figure IS^a). A similar result was ob- 
tained for all other waves. We conclude that there are no 
pathological defects in our method since we were able to 
reproduce correct output from known input. 

The second analysis provides an upper limit as to the 
contribution of leakage in this analysis. The acceptance 
is known to be significantly different than 100% and a 
function of the kinematic variables, see, for example, fig- 
ure 01 Approximating the acceptance as equal to 100% 
independent of the kinematic variables simulates a sit- 
uation where our understanding of the acceptance is as 
incorrect as possible, thus providing a limit on the effect 
of leakage. 

It was found that ignoring the acceptance does gener- 
ate spurious signals in other partial waves. Results for 
the P+ and 6*0 partial waves are shown in figure |HIb). 
Crucial to the interpretation of this spurious P4- partial 
wave is the fact that even when ignoring corrections for 
the acceptance, the A^{P+ — D+) phase difference is in- 
dependent of mass. As we will show below, the observed 



VI. EXPERIMENTAL MOMENTS 

To check the results of PWA fitting, moments of the 
spherical harmonics were calculated for each bin in m,j^o 
for each of the three t regions. These moments are la- 
belled as H{LM) and are defined as 



H{LM) = J I{n)D 



o)dn. 



Since 



H{00) = J I{fl)dn 



(5) 



(6) 



i/(00)/(47r) is just the number of events. Moments are 
directly measured observables and therefore not subject 
to ambiguities. Given a set of amplitudes, the expected 
values of the moments can be calculated. If the par- 
tial wave expansion is limited to S, P and D waves with 
\M\ < 1 there are 12 non-zero moments and the following 
relations hold: 



H{00)^\Sof + \Por + \P-r + \Do\ 



\D\ 



HilO) = —Re{SoP^} 



IP+ 



Re{PoD*} 



(7) 



+ ^Re{P^D*_} + -^Re{P+Dl} (8) 
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H{U) = l=Re{S,Pl} + -^Re{P,D*_} 

" Re{P^D*} 



Hm = l=Re{SoD;} + \ |Po|' - \ \P-f 



1 



l^ol 



1 



1 



am 



(9) 



(10) 



'^^^Re{S^D*_}-r\^-Re{P^P*J, 

+ ^Re{D^D*_} (11) 
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FIG. 9: Comparison between experimentally observed mo- 
ment distributions (filled circles) and moments calculated us- 
ing PWA solutions (open circles) as a function of 7771° effective 
mass for the moments; (a) ff(00),(b) -ff(lO), (c) H{2Q), and 
(d) ff(30). The relation between moments and PWA ampli- 
tudes and phase differences is given in the text. 



A. Fitting the moments 



^(32) = ^y|pe{P_D* } - ^y|pe{P+D;} (15) 



Re{DoD*_} 



-ff(42) 



(16) 



(17) 



(18) 



21 ' ' 21 

We note that the H{1M) and H{3AI) moments involve 
only terms linear in the P-vifave and the iJ(4Af) moments 
only involve .D-waves. The above equations also mider- 
score the inherent ambiguities in this partial wave anal- 
ysis. Solving for the amplitudes and their interferences 
in terms of the non-ambiguous observable moments does 
not lead to a unique solution. 

In figure |5| we show the comparison of moments com- 
puted directly from data with those calculated from 
the above expressions using the selected PWA solutions 
for the H{00),H{10), H{20) and iJ(30) moments. The 
agreement is excellent and similar results were obtained 
for the other moments. 



The moments as function of 7/7r° mass were fitted us- 
ing the above equations and assuming Breit-Wigner res- 
onance line shapes for the ao(980), the 02 (1320) and a 
possible P-wave resonance. These line shapes specify 
the form of the mass dependence of the amplitudes and 
phases (interference terms) in the above equations. 

A relativistic Breit-Wigner form was used: 



(m) = 



ALMBLiiqR)^) 

■ — m"^ + imrTL{m) 



(19) 



where q = X(m, 711^^,771^1) is the breakup momentum of 
mass TO. Alm is a complex, mass-independent normal- 
ization, and rrir is mass of the resonance. The energy 
dependent width is given by 



( BL{{qR?) 

\BL{{qrR?) 



(20) 



with qj. — Xinij, , rrijr , m^j) and the barrier factors Bi^{x) 
gi ven by Ba{x) = 1 , Pi (a;) = ^/x/{\ + x), B2{x) = 
\/ x'^ / {{x — + 9a;), where the R parameterizes the 
range of the interaction. In the cross section, phase space 
introduces a factor of q which multiplies the intensity 
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0.987 ±0.002 


0.980 ±0.005 


1.005 ±0.01 


Mao -' ) 


n Qfi7 -1-0 009 


QQl -1-0 007 


QQO -i-0 OOfi 

u.yyu ztu.uuo 




0.094 ±0.011 


0.246 ±0.031 


0.37 ±0.050 


Lao(ll) 


OQ/1 -1-0 01 


91 -1-0 0/1 


Q1 -1-0 0/1 




1.318 ±0.003 


1.329 ±0.002 


1.332 ±0.003 




1.318 ±0.003 


1.329 ±0.002 


1.335 ±0.003 




0.141 ±0.0093 


0.155 ±0.0067 


0.172 ±0.0092 




0.141 ±0.009 


0.147 ±0.007 


0.174 ±0.009 


Mx{I) 


1.386 ±0.032 


9.99 ±8.77 


3.47 ±0.70 


Mx(n) 


1.37 


1.37 


1.37 


rx(/) 


0.363 ±0.081 


0.11 ±0.31 


0.10 ±0.27 


rx(7/) 


0.37 


0.37 


0.37 


xm 


1.55 


1.40 


1.43 




1.54 


1.70 


2.00 



H(00) 



TABLE I: The mass and width for the ao(980), 02(1320) and 
a P-wave resonance obtained from a simultaneous fit to the 
distributions in each of 12 moments. Details are given in the 
text. For fit (I) the P-wave mass and width are allowed to 
float. For fit (II) the P-wave parameters are fixed at the values 
reported in reference fTH . The per degree of freedom, 
is also given for the two fits. 



function given in equation Q . Finally the S-wave contri- 
bution was modified by adding an incoherent background 
parametrized as a linear function in m^^o . 

There are 21 parameters allowed to vary in the simul- 
taneous fits to the distributions of the 12 moments as 
a function of rjir^ mass. The S, P and _D-wave Brcit- 
Wigners are each characterized by a mass and width for 
a total of six parameters. Each of the seven waves is char- 
acterized by a complex number but the overall phase for 
each reflectivity is fixed accounting for 12 more parame- 
ters. Finally the slope and intercept of the background 
and the R-parameter account for the final three param- 
eters. The values of these parameters resulting from the 
fits are shown in table We have also attempted to in- 
clude a second S'o-wave resonance to account for a possi- 
ble ao(1450) in the data. The moment fits were, however, 
not able to constrain the ao(1450) resonance parameters. 
We therefore removed the ao(1450) from the moments 
fits. This has no effect on the results and interpretation 
of the P-wave fits. 

The outcome of the fitting is shown for all moments 
for the three ranges of t described above in figures [TUI 
through 1211 The moments shown are corrected for ac- 
ceptance using the procedure described in reference |25l| . 
The solid curves are the result of fits in which all of the 21 
parameters described above were allowed to float. The 
dashed curves are the results of fits where the P-wave 
Breit-Wigner mass was fixed at 1.37 GeV/c'^ and the the 
width was fixed at 0.37 GeV/c^. These are the param- 
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FIG. 10: The observed H{QQ) moment as a function of 
j^tt" effective mass for different ranges of momentum-transfer- 
squared, \t^-^^^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 
H{m) = |5'o|'±|Po|'±|P-|'±|I>o|^±|£'- 



\^ + \P+\^ + \D^ 



eters reported for the P-wave exotic in the TyTr channel 
in reference 0|. 

In the low-i region, the Sq wave is prominent at m^T^o ~ 
1 GeV/c^ and agrees well with the expected contribu- 
tion to the r^TT^ spectrum from the decay of the ao(980) 
meson. The S'o-wave is produced without flipping helic- 
ities at the meson vertex and therefore it is dominant 
at low-t and decreases with increasing \t\. The structure 
at TOjj^o ^ 1 GeV/c^ in the ^^(10) and H{11) moments 
(figures ^2 and El respectively) reflects the interference 
between the S and P waves and it is noticeable only at 
low-i where the S'-wave production is the strongest. The 
wave dominates the spectrum at m,j^o ~ 1.3 GeV/c^ 
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FIG. 11: The observed H{10) moment as a function of 
777r° effective mass for different ranges of momentum-transfer- 
squared, \t^-_jj^\: (a) low \t\, (b) medium and (c) high 
The curves are results of fits described in the text. 
H{W) = -^Re{SoPS} + :^Re{PoD'o} + ^^Re{P^Dl} + 
^^Re{P+DX} 



FIG. 12: The observed H{11) moment as a function of 
■q-K^ effective mass for different ranges of momentum-transfer- 
squared, \t^-^^^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 
H{ll) = ^i?e{SoP*} + -^Re{P^Dl} - -^Re{P.Dl} 



and since it is produced with a meson hclicity flip it is 
suppressed at low-i and it is the largest in the medium-i 
region. In all three i-regions, however it is well described 
by a Breit-Wigner shape corresponding to the 02(1320) 
meson. As the D+ production increases (to medium and 
large i-regions), the structure in the moment H{\Q) (fig- 
ure 111(1 moves to the 02 region and a similar structure 
in the moment i/(30) (figure I16|l appears. Similarly, 
the interference of the negative reflectivity P- waves with 
the negative reflectivity D-waves is reflected in the mo- 
ments H{11) and Hiil) fflgures IT^ and [T7I respectivelv) 
at m^^o ^ 1.3 GeV/c^. 



The fit corresponding to the fixed P-wave parameters 
(dashed curve) is significantly worse in H{\Q) (figure ITT)) 
- a moment linear in the P-wave and somewhat worse 
or no better in the other moments that are linear in 
the P-wave - H{ll) (figure IT^ . i7(30) (figure USJ and 
i/(31) (figure [T7|l . In addition to the _ff (00) moment dis- 
tribution, the other moments - i/(20), H{22), i/(40) and 
H{A2) - are primarily sensitive to the presence of the 02 
meson, and therefore are not very sensitive to treatment 
of P- waves in the fit. The fit with the P-wave param- 
eters fixed gives a significantly worse per degree of 
freedom {xt from 1.70 to 2.00) as compared to the case 
when the P-wave parameters are allowed to fioat (x^ 
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FIG. 13: The observed H{2Q) moment as a function of 
r/Tr" effective mass for different ranges of momentum-transfer- 
squared, \t^-_^^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 
H{20) = ^i?e{5oA1}+| |Po|'-||P-p-i IP+P + I 
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D4 



FIG. 14: The observed H{21) moment as a function of 
■q-K^ effective mass for different ranges of momentum-transfer- 
squared, \t^-^^^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 



from 1.37 to 1.43). In the later case, however, it was 
found that no single set of consistent parameters could 
be found for the P waves across the three i-ranges. In 
the \aw-t range the P-wave with mass of 1.39 GeV/c^ 
and width of 0.36 GeV/c^ was found, which is the only 
solution consistent with what was reported for the ryTr^ 
channel However, even in this t-range other solu- 

tions could also be found with masses ranging from 1 to 
5 GeV/c^. For example, a solution for a P-wave with a 
mass of 0.91 GeV/c^ and a width of 1.05 GeV/c^ had 
a, xi = 1.57. This is because the interference with the 
ZJ-waves is weak at low-i and the S'-wave is too narrow 



to constrain the P away from the m^TrO ~ 1-3 GeV/c^ 
region. The medium-i to large-t range should be the 
best for constraining the P waves. Here the D waves 
are dominant and some S'-wave is still present. In this 
t-range, however, acceptable solutions (with comparable 
xi) different from the one listed in tabled] could not be 
found. 

The non-vanishing moments H{1Q) f figure [TT|l . -ff(ll) 
(figure [121, H{m) (figure CHll, 77(31) (figure [T2|) and 
H{32) (figure [TH|l clearly indicate a presence of a P-wave 
in the jyTr*^ system. The large values of mass and width 
found indicate that the Breit-Wigner parameterization 
for these waves is most likely inadequate. 
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FIG. 15: The observed H{22) moment as a function of 
r/Tr" effective mass for different ranges of momentum-transfer- 
squared, \t^-_^^\: (a) low (b) medium \t\, and (c) high 
The curves are results of fits described in the text. 



FIG. 16: The observed H(30) moment as a function of 
■q-K^ effective mass for different ranges of momentum-transfer- 
squared, \t^-^^^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 
Hm = |y|i?e{PoI?a-7^^?e{P_i31}-^i?e{P+7?;} 



VII. MASS AND t DEPENDENCE OF PARTIAL 
WAVES 

A. PWA solutions and fits to the moments 

In figures through [21 various partial waves and 
phase differences are shown as a function of rji:^ mass for 
the low-|i|, medium- |i I and high-|i| ranges. In particular 
we show the mass dependence of the Sq and Dq ampli- 
tudes and the A<I>(S'o — Dq) phase difference. Also shown 
is the mass dependence of the P-y. and Z?+ amplitudes 
and the A<I>(P-|_ — £)+) phase difference. The filled data 
points indicate the selected (physical) solutions. The ad- 



ditional points indicated by open circles and dashed error 
bars are those ambiguous solutions that satisfy conditions 
described below. 

The curves are the result of the fits to the moments 
where the P-wave Breit-Wigner parameters are allowed 
to float. We refer the reader to similar plots shown in 
figurc[3for the and D+ amplitudes and the A<I>(S'o — 
Dq) phase difference along the the curves used to select 
the solutions. 

The results of the fits to the moments (recalling that 
the moments are unambiguous) were used to study the 
possible inclusion of other ambiguous solutions not se- 
lected as the physical solutions. We compute a that 
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FIG. 17: The observed ^(31) moment as a function of 
r/Tr" effective mass for different ranges of momentum-transfer- 
squared, \t^-_jj^\: (a) low \t\, (b) medium and (c) high 
The curves are results of fits described in the text. 



FIG. 18: The observed H{32) moment as a function of 
rjiv'^ effective mass for different ranges of momentum-transfer- 
squared, \t^-^ji^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 



measures the deviation of other ambiguous solutions from 
the partial wave solutions obtained from the moments 
fits. This requirement was imposed separately for the 
^0 and D_|_ -waves and the A<l'(5o — Dq) phase differ- 
ence. In some cases several ambiguous solutions were 
found to have similar to the selected solutions and 
these are included in figures figures through 1241 as the 
open circles with dashed error bars. The mass-dependent 
fits involving only the P-f and amplitudes and the 
A$(P+ — - as described below - phase difference 
were insensitive to the inclusion of these solutions. 

The fits shown in figures through |^ describe the 
Dq and in all three \t\ regions. The values obtained 



for the mass and width of the a2(1320) as shown in table 
^are consistent with nominal values for this well under- 
stood state. The Sq is well described in the low-|i| region 
and the mass and width shown in table for the ao(980) 
agree with published values. As |i| increases the produc- 
tion of the flo (980) decreases and the description of the 
Sq wave in figure !^ for the high-|t| range is problematic. 
In this region the 5*0 wave includes mostly structureless 
background naturally subsumed into the Sq wave. Re- 
quiring that the wave be described by a Breit-wigner 
wave yields a width that is wide and driven by the phase 
motion of the a2(1320) via the A<I>(S'o — Dq) phase dif- 
ference. As noted above, the moment fits did not con- 
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FIG. 19: The observed _H'(40) moment as a function of 
r/Tr" effective mass for different ranges of momentum-transfer- 
squared, \t^-_jj^\: (a) low \t\, (b) medium and (c) high 
The curves are results of fits described in the text. 



FIG. 20: The observed -^(41) moment as a function of 
rjTv'^ effective mass for different ranges of momentum-transfer- 
squared, \t^-^^^\: (a) low (b) medium and (c) high 
The curves are results of fits described in the text. 
H{41) = lyqRe{DoD*_} 



strain the parameters of a possible second scalar reso- 
nance, the ao(1450). In the mass-dependent amplitude 
and phase difference fits, the inclusion of a second scalar 
resonance In the low-i region yields M = 1.34±0.06 GeV 
and r = 0.41 ± 0.09 GeV for the mass and width of the 
second qq resonance, while in the medium-i region we 
obtain M = 1.44 ± 0.03 GeV and T ^ 0.25 ± 0.06 GeV, 
respectively. No stable fit was found for high-t data. 

The moments fits poorly describe the P+ waves in all 
three \t\ regions and the P+ amplitude peaks in the vicin- 
ity of the dominant wave strongly suggesting leakage. 
We also note that the P+ wave has an additional lower 
mass peak in the vicinity of the ao(980) in the low-|t| 



region (see figures |221 and |23l where production of the 
ao(980) is strong - again suggesting leakage. Our Monte 
Carlo studies indicate that leakage of dominant waves can 
lead to false peaks in other amplitudes but it is unlikely to 
cause false phase motion. Indeed this P+ wave enhance- 
ment at about the mass of the ao(980) follows closely 
the behavior of the S wave in its dependence on both 
mass and momentum transfer squared. Furthermore as 
shown in figure [221 and especially in figure the addi- 
tional solutions allowed by the selection criteria (those 
indicated by open circles) do not alter the S wave shape 
but reduce the amplitude of the P+ wave in the region 
of 1 GeV/c^. Furthermore, the fits to the moments 
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FIG. 21: The observed H{A2) moment as a function of 
r/Tr" effective mass for different ranges of momentum-transfcr- 
squared, \t^-_^^\: (a) low \t\, (b) medium and (c) high 
The curves are results of fits described in the text. 
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do not reveal any P w&yb signal at ~ 1 GeV/c^ and due 
to angular momentum barrier effects it is expected that 
at low breakup momentum the S wave should dominate 
over the P wave. 



> 

a; 
O 



3000 


_' 




2500 






2000 






1500 






1000 






500 














1500 
1000 
500 



0.8 1 1.2 1.4 1.6 1.8 2 



T]X effective mass (GeV/c ) 




0.8 1 1.2 1.4 1.6 1.8 2 

2 

T]X effective mass (GeV/c ) 



FIG. 22: PWA solutions for the: (left) - (top) So (middle) Dq 
waves and (bottom) A$(S'o — i3o) phase difference and for the 
(right) - (top) P+ (middle) D+ waves and (bottom) A$(P+ — 
D+) phase difference for the low-|t| region. The fits are result 
of the fit to the moments. The filled diamonds correspond to 
the selected solution. The open circles are additional solutions 
allowed by the criteria described in the text. 



Wigner distributions corresponding to the resulting fits 
are shown in figure [23 for the Z?+, ZJq, D- waves and 
the two Dm=2 waves of opposite reflectivities (the two 
M — 2 curves are indistinguishable in the plot). This 
supports our original assumption. 

The systematics of varying all possible fit assumptions 
was studied for the Z3+-wave for the high-|i| region and 
the results are shown in the plot of figure[5B| We chose to 
study this high-|f| region since the description of the 
wave amplitude is poorest for this region ((see figure 121 • 
The curve is the result of a of the overall fit to the mo- 
ments. The shaded region corresponds to the variations 
in the fit by letting the Breit- Wigner parameters for all 
the other waves vary in all possible combinations - that 
is, being fixed at nominal values or allowed to float. The 
ZJ-i- wave is within the systematic errors. 



B. Systematics and the fitting procedure 

The assumption was made earlier that waves with 
|M| > 1 were not important. In order to test this as- 
sumption the additional \M\ — 2 D waves, one for each 
reflectivity, were included in the fits for all the moments 
including the iJ(43) and iJ(44) moments. The Breit- 



C. Mass-dependent fits 

The dependence of the measured and P_|_ ampli- 
tudes and the A$(P_|_— £>+) phase difference as a function 
of rjTT^ mass was fitted to two interfering Breit- Wigner 
line shapes following the method of reference ^Jl- Fits 
carried out separately for the three ranges in t are pre- 
sented in the plots of figure and the fit parameters 
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I ' I I 1. u 1^ — I — I I ^ — I 

0.8 1 1.2 1.4 1.6 1.8 2 0.8 1 1.2 1.4 1.6 1.8 2 



T]X effective mass (GeV/c ) T]X effective mass (GeV/c ) 

FIG. 23: PWA solutions for the: (left) - (top) So (middle) 
Do waves and (bottom) A$(S'o — Do) phase difference and 
for the (right) - (top) P+ (middle) D+ waves and (bottom) 
A"I>(P+ — D+) phase difference for the medium-|t| region. The 
fits are result of the fit to the moments. The filled diamonds 
correspond to the selected solution. The open circles are ad- 
ditional solutions allowed by the criteria described in the text. 




0.8 1 1.2 1.4 1.6 1.8 2 0.8 1 1.2 1.4 1.6 1.8 2 



T]j: effective mass (GeV/c ) T]j: effective mass (GeV/c ) 

FIG. 24: PWA solutions for the: (left) - (top) So (middle) Do 
waves and (bottom) A$(S'o — -Do) phase difference and for the 
(right) - (top) P+ (middle) D+ waves and (bottom) A<1>(P+ — 
D+) phase dilference for the high-|t| region. The fits are result 
of the fit to the moments. The filled diamonds correspond to 
the selected solution. The open circles are additional solutions 
allowed by the criteria described in the text. 



are listed in table ^1 In figure |2H1 the P-wave param- 
eters obtained from the mass dependent fits are com- 
pared to those reported in references i2, 13]. The 
P-wave mass obtained for the high-i range is consistent 
with that reported in the r]TT~ channel but the width ob- 
served here is significantly higher. The values obtained 
for the mass in the low-i and high-i ranges are lower than 
for the r^TT" channel. The conclusion is that no consistent 
set of resonant P-wave parameters can describe the data 
while the resonant parameters obtained for the a2(1320) 
are consistent for various ranges in TyTr" mass and \t\ and 
consistent with well-established parameters from earlier 
experiments. 



D. Description of the t dependence of the D- waves 

Figure I^Hl shows the dependence of the , Dq and 
I^.-waves as a function of t^-^^^. Since the 02(1320) 
is dominant in these waves and the production of these 
waves is associated with exchanges well known from 
Regge phenomenology, a comparison of data with this 
known phenomenology is a check on the methodology. 
Indeed the t dependence is well described. 

The wave is produced via a natural t-channel 
charge exchange which is dominated by the p trajectory 
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FIG. 25: Comparison of the intensity of the D waves for the 
medium-|t| region including the D waves with M—2. There 
are two M=2 waves with opposite reflectivity. 
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FIG. 26: The D+ wave as a function of 7771'' effective mass 
for the high-|t| range. The curve is the result of an overall fit 
to the moments, including a D-wave resonance, as described 
in the text. The shaded area is an estimate of the system- 
atic errors derived from a variation of all combinations of the 
Breit-Wigner parameters for the 5", P and D amplitudes - as 
described in the text. 





all |t| 


low-|t 


medium- t 


high-|i| 




1.326 


1.316 


1.329 


1.326 




±0.0023 


±0.0049 


±0.0029 


±0.0036 




0.169 


0.127 


0.154 


0.166 




±0.0069 


±0.014 


±0.0082 


±0.01 


Mx 


1.272 


1.301 


1.268 


1.356 




±0.017 


±0.014 


±0.023 


±0.021 


Tx 


0.66 


0.19 


0.67 


0.629 




±0.048 


±0.032 


±0.087 


±0.064 


xi 


3.23 


2.13 


1.51 


1.60 



TABLE II: The masses and widths obtained from a simultane- 
ous fit to D+ intensity, the P+ intensity and the A<I>(P+ — _D+) 
phase difi'erence. For these fits both waves are parametrized 
as Breit-Wigner resonances. In addition to the fits for the 
three ranges of t we also include the fit for the full t range to 
compare with the results of reference [Tl|. 
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FIG. 27: PWA solution for the (top) P+ (middle) D+ waves 
and (bottom) A$(P+ — D+) phase difference for low-|f|, 
medium-|t| and high-|f| as a function of ij-rv^'' effective mass. 
The curves are the result of a mass-dependent fit assuming 
Breit-Wigner resonance forms for the P+ and D+ waves. 
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FIG. 28: Comparison of P-wave mass and widths among var- 
ious studies. A mass dependent fit to the D+ and P+ am- 
plitudes and A"1>(P+ — D+) phase difference, assuming Breit- 
Wigner line shapes was carried out for the E852 data for tj-k" 
- this study - and Tyvr". Results are also compared to the 
Crystal Barrel results. 



and is parametrized by, [2£ 



D 



++ 



^bpt^-i-Kap{t)/2 ^ 



r(l-ap(t))sin(™p(i)/2) (21) 



(22) 




Here the and D^~^ are the nucleon helicity flip 

and non-flip amplitudes respectively, with r = 7 in equa- 
tion [23 according to reference |2^ , and 



(23) 



The parameters of the p trajectory are given by ap(t) — 
0.5 + 0.9t [GeV^], bp = 3.25 GeV'^ [s^. The flt of 
the i-distribution given by Eq. H22|) to the measured t- 
distribution involves one parameter - the overall normal- 
ization and is shown in flgure 123 Some discrepancy is 
seen at large- t, in particular the characteristic dip of the 
p exchange at i = is not seen. This indicates that 
there is some absorption required for the 02 . 

The D'^ wave is produced via unnatural exchange 
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5000 




t' (GeV/c) 

FIG. 29: Dependence of D+ (circles), Do (triangles) and D- 
(diamonds) waves on \t'\ or \t^-_^^^\ — \tmin\- The curves are 
the results of fits to forms for these distributions expected 
from Regge phenomenology. 



and is expected to be dominated by the &i in the nu- 
cleon helicity flip amplitude and the tensor, Z-trajectory, 
(J^*-^/ = 2 1) in the non-flip amplitude. These can 
parametrized by, 



D+-=^ 



1 



t 



-iTTOb-^ (t)/2 



X 



D++ 



r{-ab, (i))sin(^ahi W/2)(24) 
fcztg-^^az(t)/2p(^ _ cos(7raz(t)/2). (25) 



with the parameters of the two trajectories given by 
ab, = -0.37 + 0.9f, [GeV^]^z = 0.9t, [GeV^], and 
^bz^b^^ 4.51 GeV-2 ^ . 
The Z3_-wave is produced by unnatural parity ex- 
change and its t dependence is given by D^{t) — 
^y—t'/4:m^Do{t). The fit to the measured Dq wave in- 
tensity has two parameters - the normalization of the 
nucleon helicity flip and non-flip amplitudes - and is also 
shown in figure 1291 



E. Ratio of D-waves 

In a paper discussing the PWA of the rjir system pro- 
duced in peripheral np interactions l3ll|, Sadovsky uses 
a criterion for selecting the physical solution among the 
ambiguous solutions. The test involves a measurement of 
the ratio rg — {Dq — D_)/D^. According to Regge phe- 
nomenology this ratio should fall like a/p-jr where a is de- 
termined from other experiments. At Pti=18 GeV/c this 
ratio is expected to be rg = 0.80. Integrating the distri- 
butions in figure 1^ over all t we obtain rs = 0.72 ±0.12. 



VIII. CONCLUSIONS 

A partial wave and moment (averages of spherical har- 
monics) analysis of data collected in experiment E852 
from the reaction n^p rjir^n (where — > 77) as a 
function of rjir'^ effective mass and momentuni-transfer- 
squared was performed. The presence of two well- 
established ao(980) and a2(1320) resonances allows for 
the imposition of criteria for selecting the physical so- 
lution among the mathematically ambiguous solutions. 
The physical solution is consistent with the presence of 
two interfering Breit-Wigner resonances for these states. 
The distribution in 12 acceptance-corrected moments, 
H{L,M), calculated directly from the data are consis- 
tent with moments calculated from the PWA solutions. 

Every term in the moments with odd L (5 of the 12) 
is linear in the P wave and the observed moments dis- 
tributions are non-zero clearing indicating that the data 
demand the presence of a P wave. To allow for the inclu- 
sion of a possible resonance with exotic quantum num- 
bers, the moments are thus fit with the assumption of 
three interfering Breit-Wigner resonances for the S'-wave 
and 13- wave and P-wave. 

The data suggest weak production of the ao(1450). 
The ao(1450) was established by the Crystal Barrel col- 
laboration in pp — !■ ir^TT^r] and pp Tr^r]r] |28j | with pro- 
duction ratios for ao(1450) : ao(980) : a2(1320) equal to 
4 : 6 : 33 and 11 : 14 : 25 respectively. These are quali- 
tatively consistent with our results. The dominant 13+ , 
Do and D.-waves are consistent with production of the 
02(1320) and fits to the their line shapes are consistent 
with the resonance parameters of this well-established 
tensor resonance. Moreover, the observed t-dependence 
of the individual D waves are well described by Regge 
phenomenology as is their ratio. Finally the qq and a2 in 
the data are well described as Breit-Wigner resonances 
produced via a t-channel exchange mechanisms regard- 
less the treatment of the P-wave, and/or incoherent back- 
ground. 

While the Breit-Wigner parameterization of the S 
wave and D wave is physically well motivated and con- 
firmed by the data, the resonance representation of the 
P-wave is problematic. The data unambiguously indi- 
cated that the phase of the P-wave increases as the func- 
tion of the rjTr'^ mass which corresponds to an attractive 
interaction in the rjir'^ P-wave channel. However, the 
data do not require the phase of the P-wave to increase 
over 90° nor is there compelling evidence for the inten- 
sity of the P-wave to have a resonance shape in the mass 
range studied, i.e from threshold to 1.8 GeV. 

Mass-dependent fits, restricted to the P+ and 13+ am- 
plitudes and A<I>(P+ — £)+) phase difference, following the 
method of reference jlJI reporting an exotic J^'-' = 1 ^ 
meson decaying into T77r~, yield results inconsistent with 
those of reference [ll| and among the three t ranges in 
this analysis. The discrepancy in the resonance parame- 
ters found by fitting the natural waves alone is, however, 
smaller than the one found by fitting the full data sample 
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(moments). 

However, by restricting the study to only P+, D+ 
waves and the phase, Ai>(P+ — £*+), important addi- 
tional information from the remaining data is ignored. 
For example, an exotic P wave resonance should have 
the same resonance parameters in all three, P-,Po and 
P+ -waves. These wave represent different spin projec- 
tions of the resonance and the data should discriminate 
between their production strengths. 

Furthermore, we have found that the intensity of the 
weak waves, i.e. P and S'-wave above the 0,0 region are 
more affected by leakage from the strong, D waves due 
to the 02 than are the phase differences. Thus, since the 
mass-dependent fit to the natural waves alone, is strongly 
constrained by the intensity distribution of the P-wave, 



it tends to follow the 02 line shape, and force the fit to 
mimic a P-wave resonance with mass near the 02 mass. 

In a separate paper we will present an alternative de- 
scription of the mass dependent P-wave amplitude and 
phase that does not require the existence of an exotic 
meson but is consistent with rjir re-scattering. 
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